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Abstract 



We use (nonconservative) dynamical semigroups to investigate the decay 
law of a quantum unstable system weakly coupled with a large environ- 
ment. We find that the deviations from the classical exponential law are 
small and can be safely ignored in any actual experiment. 
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The decay law of an unstable microscopic system can be rather well described by an 
exponential function; this result can be easily justified on the basis of classical probabilistic 
considerations. However, microscopic systems should be described by quantum mechanics 
and it is well known that in quantum theory the exponential decay law can not be valid 
for all times; in particular, it surely fails for very short and very long times. Indeed, rather 
general considerations assure that the quantum decay law can be described by a three-step 
function: a Gaussian law at short times, the classical exponential law at intermediate times 
and finally a power law at longer times (e.g. see [1-3] and references therein). 

Although various experiments have been devised in order to obtain evidence for dis- 
crepancies with the exponential law, none have been so far actually detected. [4] In this 
respect, a system that has attracted a lot of interest both theoretically and experimentally 
is the neutral kaon system. The K°-K° system has proven to be one of the most fruitful 
systems for testing fundamental symmetries, like CP and CPT. In parametrizing viola- 
tions of these symmetries, one usually takes for granted the exponential decay law and 
uses an effective theory to describe the kaon system. [5] 

In the following we shall examine to what extent the exponential decay law can be 
trusted in actual experiments. More precisely, we shall study the deviations from the 
exponential law for small and large times, taking also into account possible effects due to 
incoherent interactions with the environment. 

The general idea that is at the basis of our considerations is that unstable systems 
can be viewed as specific examples of open quantum systems. These systems can be mod- 
eled in general as being small subsystems in weak interaction with large environments. 
Although the global time evolution of the closed compound system is described by an 
unitary transformation, the reduced dynamics of the subsystem, obtained by the elimi- 
nation of the environment degrees of freedom, usually develops some sort of dissipation 
and irreversibility. Under mild assumptions, the reduced evolution is realized by one- 
parameter (=time) maps acting on the states of the system, conveniently described by 
density matrices, with forward in time composition (semigroup property) and the addi- 
tional characteristic of being completely positive. This set of transformations forms a 
so-called dynamical semigroup. [6-10] 

This rather universal and general formalism has been recently adopted to treat effec- 
tive dynamics for the kaon system that transform pure states into mixed ones. [11-16] The 
physical motivations behind such an approach are based on quantum gravity, that predicts 
loss of quantum coherence at the Planck's length due to fluctuations of the gravitational 
field. [17] These generalized time-evolutions lead to CP and CPT violating effects that 
could be in the reach of the next generation of neutral kaon experiments. [15, 16] 

In these treatments, the environment was assumed not to contribute to the decay 
of the kaons, that was effectively described by the standard exponential law. However, 
possible effects of the environment on the decay process itself are surely conceivable. As 
we shall see, they can be studied using again the formalism of dynamical semigroups. It 
will turn out that these environment effects do not modify the short-time behaviour of the 
decay law, but only the exponential and power law regimes. 

However, in the realistic hypothesis of a weak coupling between unstable system and 
environment, these modifications are tiny or occur for too large times for any practical con- 
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siderations. Therefore, the exponential decay law captures rather well the essential features 
of the decay process of an unstable system and can be certainly used with confidence in 
any actual experiment involving weak-decaying particles, like the neutral K mesons. 

For our considerations, we choose to describe the states of a quantum system evolving 
in time by means of density matrices p(t). Given an initial state p and a time-independent 
hermitian hamiltonian H, the standard time evolution is described by 

p(t) = e~ iHt p e lHt , (1) 

solution of the Liouville-von Neumann equation 

In the case of an unstable system, it is custom to split H as 

H = HW + HW, (3) 

where is the unperturbed hamiltonian, while H^- 1 ' is the interaction hamiltonian that 
drives the decay process; the system would be stable if = 0. In the case of the neutral 
X-mesons, 

#(0) 

can be identified with the hamiltonian of the strong interactions and > 
with that of the weak interactions (possible "superweak" mixing terms should be included 
m HW). Further, we shall call P u the projector operator on the subspace of the undecayed 
states. We shall also use the orthogonal projector = 1 — P u ; it describes the transition 
to the space of the decayed states. 

Assuming that at the beginning our unstable system is in the undecayed state p = 
P u P Pu, the probability of finding it undecayed at time t is given by 

V(t) = Tr[p(t) P u ] . (4) 

In the case of the decay of a single particle, with evolution as in (1), the properties of 
V(t) have been widely studied in the literature. Here we shall extend those treatments by 
considering also the interaction of the decay system (not necessarily one-dimensional) with 
the environment. Together with the effects of the hamiltonian also this interaction 

could in principle contribute to the decay process. 

As mentioned in the introductory remarks, in order to take into account these extra 
effects we shall treat the unstable system as an open quantum system S. [6-10] As for 
any open system, S can be thought of as interacting with a suitable environment E, so 
that the global system S + E is closed. This evolves in time according to a group of 
unitary operators as in (1), governed by a total hamiltonian Hs+e, which is the sum of 
the hamiltonian Hs of the subsystem, of the hamiltonian He of the environment and of the 
interaction hamiltonian Hs e between them. A reduced dynamics for the subsystem S can 
be consistently obtained when S and E are assumed to be uncorrelated at the moment of 
the formation of the unstable system and therefore the state of the total system is simply: 

PS+E = P® PE- 
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(2) 



In such cases, by tracing over the the environment degrees of freedom one gets linear, 
completely positive maps p %[p] on the states of the subsystem, where 



lt[p] = Tte e 



iHs + Et (p(g) p E )e 



,iH s +Et 



(5) 



These transformations do not have any simple composition law and they usually contain 
memory effects. However, when the interaction with the environment is weak, they reduce 
to completely positive, dissipative dynamical maps 7$ : p 1— > p(t), which obey a semigroup 
composition law, r y tl 07^ = 7t 1 +t 2 , t±, t<i > 0, and satisfy the condition of entropy increase: 
dS/dt > 0, S(t) = — Tr[p(t) lnp(f)]. Moreover, these semigroups are generated by equations 
of a very specific type that can be explicitly given: [6] 



where the set of (bounded) operators R and A k are such that: R + R j + J2k A k A k < 0. 
They are probability preserving or not depending on whether this combination is identically 
zero or not; they are called quantum dynamical semigroups. Notice that this description of 
open systems is rather general and the effective subdynamics j t is essentially independent 
from the type of the environment. 

In the case of unstable high energy particle systems, a natural choice for E could 
be the gravitational field, whose effects are usually neglected in the theory of elementary 
particles because of the smallness of its coupling. Its quantum fluctuations at Planck's 
length could nevertheless act as a weak coupled environment producing detectable effects 
in elementary particle interactions. [17, 11] It is also worth mentioning that the effects of 
such fluctuations, realized via the space-time foam, can be effectively described with an 
heat bath, the most natural of all environments. [18] 

In view of all above considerations, we shall now generalize the standard quantum 
mechanical evolution equations for an unstable system by adding to (2) a linear piece of the 
form given by the r.h.s. of (6). Since this additional piece should describe the contribution 
of the environment to the decay process, it must be proportional to pit) — P u p(t) P u . In 
fact, this term would project p(t) out of the space of undecayed states. Thus, we shall 
study the generalized time evolution described by the equation 



with A a positive "coupling" constant and H as in (3). We shall consider weak coupled 
environments, and therefore assume A to be much smaller than any typical energy scale 
in the hamiltonian. In the case of the K°-K° system, assuming the dissipative term 
in (7) of gravitational origin, dimensional arguments suggest A to be at most of order 
m 2 K /mp, where rriK is the kaon rest mass and mp is Planck's mass. The evolution equation 
(7) describes nonconservative dynamical maps. [6] Indeed, probability is not conserved, 
(d/dt)Tr[p(t)] = —XTr[p(t) Pd] < 0, as it should be for an unstable system; however, this 
violation is small, since A is small. 





k 



dp(t) 



i[H,p(t)]-\(p(t)-P u p(t)P u ) , 



(7) 



dt 
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Equations of the form (7) were used before in discussing one-dimensional unstable 
systems. [1] The motivation for the introduction of the nonstandard term was there at- 
tributed to the interaction of the decaying system with the measuring apparatus. While 
this point of view is surely viable, we stress that our interpretation of (7) as describing the 
evolution of an unstable system is completely phenomenological in nature. In particular, 
we do not make any specific assumption on the phenomena responsible for the appearance 
of the second term in the r.h.s. of (7). 

In order to proceed in the study of equation (7), it is convenient to introduce a vector 
notation, rewriting the density matrix p describing the state of the unstable system as the 
vector \p). To any operator acting on the state p, one can define a corresponding operator 
acting on the vector \p). In particular, one can define the projector LT U , 

n u \p) = \p uP p u ) , n 2 u = n u , (8) 

and its orthogonal complement: lid = 1 _ LI U . By introducing the Liouville operator Lh 
corresponding to the hamiltonian H, 

L H \p) = \[H,p]) , (9) 

one can rewrite (7) as a Schrodnger like equation: 

i^\p(t)) = L \p(t)) , (10) 

with 

L = L H -i\n d . (11) 

Although the statements and the conclusions obtained below using this new formalism can 
be rigorously justified, for sake of simplicity we shall keep mathematical considerations to a 
minimum. As we shall see in the following, this formalism results particularly appropriate 
in the study of unstable system for which the space of the undecayed states is not one- 
dimensional, as in the case of the neutral kaons. 

We are ultimately interested in describing the properties of the probability V(t) in 
(4) . Therefore, one should concentrate on the study of the time evolution of the projected 
vector: 

\ P (t)) u = n u \ P {t)) . (12) 

Using the Laplace transformed vector \p(s)), the corresponding Lippmann-Schwinger equa- 
tion reads: 



s + i I L uu + L u d - ■ — . . Ldu 

is - L dd + i\n d 



\P{s))u = W))u , (13) 



where s is the Laplace variable and 

n u L H n u , L ud = n u L H n d , L du = LidL H n u , L d d = n d L H nd . (14) 



J uu 
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The action of these operators on any state \p) is well-defined. For instance, using the 
shorthand notation Oij = Pi O Pj, i,j = u, d, with O a generic operator, one finds: 



L u d \P) = \H u dPdu — PudHdu) 
Ldu | p) = | Hdu Puu — Puu H u d ) 



(15a) 
(156) 



According to our hypothesis, the decay process is driven by the interaction hamiltonian 
and the coupling to the environment. In particular, the decaying system would be 
actually stable for = and A = 0. This immediately implies that H ud = H ud \ 

Hdu = H^J and H uu = H u °u . Furthermore, as is supposed to be small (in the case 
of the neutral kaons, H^- 1 ' is indeed the weak hamiltonian), we shall study approximated 
solutions of (13), taking into account only terms up to second order in H^; this is also the 
approximation that is usually adopted in the description of decaying particles in standard 
quantum mechanics. 

At this point, in order to simplify the formulas, we also assume that the space of the 
undecayed states is degenerate in energy: H uu = H u °u = EqP u . In the case of the neutral 
kaons, this is not really a restriction since K and K° have the same rest mass. Then, by 
acting on a generic state \p) u , one can prove that: 



Luu L t 



1 



J du 



-W{s) 



+ o 



i(s + \)n d - L dd 
where W(s) is the effective non-hermitian "hamiltonian" 

W( s ) - H W * 

W \ S >- n ud .. TT ( ) n du > 

and ^vk(s) is the corresponding generalized Liouville operator 

C W (s)\p)= W{s)p-pW\s) 



((# (1) ) 3 ) , 



(16) 



(17) 



(18) 



This is not surprising since, in view of (15), the generalized operators L u d arid Ldu them- 
selves are of order Thus, up to second order terms in equation (13) becomes 



(19) 



s + iC w{s) ] \p{s)) u = |p(0)), 
Using the inverse Laplace transform, one can then write 



/•c+ioo 


1 


/ ds e st 


_s + i£w(s) . 


I c—ioo 



\P(0))U , 



(20) 



where c must be chosen so that the integration path in the complex s-plane lies to the 
right of all singularities of the generalized operator [s + iCw{ s )]~ l - 
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Before proceeding further, let us first deduce from (20) the standard exponential decay- 
law. The singularities of the integrand in (20) are related to those of the Liouville operator 
£w(s)i an d therefore to those of W(s) in (17). Indeed, using the definition (18), one can 
deduce that the spectrum of £w(s) coincides with the difference of the spectra of W(s) and 
W* (s). The operator W(s) is analytic in the complex s-plane except when the denominator 
in (17) vanishes. Since typically the spectrum of H^J is continuum, the singularities of 
W(s) lay on the imaginary axis, where usually there is a cut. Therefore, one can conclude 
that also [s + iC\y^]~ 1 has generally a cut on the imaginary axis of the s-plane. One can 
then move the integration path in (20) almost to coincide with this axis. 

Further, one notice that when H^- 1 ' = and A = the generalized Liouville operator 
£\y(s) becomes the null operator, so that only a pole at the origin occurs in the integrand 
of (20). In this case, we obtain: \p(t)) u — |p(0)) u , as it should be, since now the system is 
stable. In presence of interactions, this pole moves into the second sheet, but remains close 
to the imaginary axis for small and A. [3] One then expects that this pole continues 
to give the main contribution to the integral. This is the so called Weisskopf-Wigner 
approximation, that gives rise to the exponential decay law for all times. Indeed, within 
this approximation and using (18), the integral in (20) gives 

\p(t)) u = e~ itc »w \ p (0)) u , (21) 

or equivalent ly 

Puuit) = e~ iHwt p uu (0) e m w f , (22) 
where the effective hamiltonian H\y takes the form: 

Hw = H<$ + H% 1 iff • (23) 

This hamiltonian is not hermitian and can be written as Hyy = M — iT/2, with M and r 
hermitian and positive. Indeed, by putting the system in a finite box so that the spectrum 
of H^J becomes discrete, one explicitly finds: 

iMU = E 5 a p + J2(*\H {1) \k) {E ^ Ek yl A 2 (k\H (1) \P) , (24a) 

[r] a0 = 2j2(*\H^\k) _ ^ k)2 - A2 (k\H^\[3) , (246) 

where we have used greek (latin) indices to label undecayed (decay-products) states. The 
entries of the matrix 

[a(E)U= (a\H {1) \k) (k\H^\(3) , (25) 

E k <E 

are usually found to be piecewise differentiable in the variable E when the volume of the 
box becomes infinite. [2] In this limit, by setting 

_ d[a(E)] a p 

K-kJW - — ^ ' ( 2t V 
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the sums in (24) can be substituted by integrals: 



[M} aP = E 8 a p 



dE [u{E)U 



(E - E) 



(E - E) 2 + X 2 



/>oo 

[ru = 2 / dE [u{E)] aP -—- 



A 



E) 2 + X 2 ' 



(27a) 
(276) 



where E m is the lowest eigenvalue of H^ d . For A small, using 

A =7t5(x)+Ap4+°( A2 ) ' 



a;- 



+ A 2 



ar 



with P indicating principal value, the equations (27) finally become: 

[M] Q/3 = [v{E )] af3 + Att {u'(E )] af3 , 

[r] a p = 2tt [u(E )U + 2A (<W - W(Eo)U) , 

where the dash signifies derivative with respect to E and 

[//(£o)W = E S aP + / dE [u(E)] al3 P- 

JE m E - E 



(28a) 
(286) 



(29a) 
(296) 



(30) 



When A = 0, M and r reduce to their standard Wiesskopf-Wigner expressions. The effect 
of the environment is to modify these expressions by adding terms that are in principle 
calculable using field theory techniques. [19-21] The actual evaluation of the various terms 
in (29) requires the adoption of a specific microscopic model for the interacting hamiltonian 
and is certainly beyond the purpose of the present work. 

As stressed at the beginning, the exponential decay law can not hold for all times. 
We shall now go back to the equation (19), which is exact up to second order terms in 
H^ l \ and try to evaluate the corrections to the Weisskopf-Wigner approximation. To this 
purpose, let us add and subtract to the l.h.s. of (19) the term i£# w .[22] Within our 
approximation, one can then write: 



s + iC 



W(s) 



(31) 



By acting on a generic state \p) u , using simple algebra one can prove that: 

-^2 (f-W{s) - C-Hw) \p)u = | V(s) p uu - p uu V^s) ) , 



where 



V(s) = H {1) 



1 



(Eo-HZ'+iXr 



1 



i(s + \) + E -H i 



(0) is 
dd 



n du 



(32) 



(33) 



Inserting these results in (20) and performing the s-integration, one obtains an effective 
evolution for the projected density matrix p uu of the form 



Puu{t)=U{t) Puu (0) U\t) 



where 



U(t) 



-iH w t _|_ 



(1) 
ud 



i(H ( fJ -i\)t _ -iE t 



H 



(i) 

du 



(34) 



(35) 



This evolution satisfies the correct boundary condition, U(0) = P u ; it contains contribu- 
tions up to second order in and to all orders in A. It is dominated by the exponential 
Weisskopf-Wigner term; however, the additional pieces are relevant for small and large 
times. 

By expanding (35) for small t, one gets: 



U(t) ~ (1 - itE )P u - - (e%p u + h£1 h£>) 

Inserting this result in (34), the probability (4) takes the form: 

V(t) ~ 1 - (t/r G ) 2 ee 1 - t 2 Tt\h£1 H$ Puu (0) 



(36) 



(37) 



and therefore has a Gaussian behaviour; the Gaussian width Tq 2 , that can be rewritten 
as Tr[H u d Hdu Puu(ty], gives the spread in energy (AE) 2 of the initial state p uu (0). Fur- 
thermore, V(t) is independent from A, so that the small-time decay law is unaffected by 
the interaction with the environment. This result is physically understandable: the time is 
too short to allow the environment to play a role in the decay, which, in this early stages, 
is totally driven by H^ 1 ' . This is also in agreement with our starting assumption that the 
environment does not disturb the preparation of the decaying system. 

The times for which the Gaussian behaviour of V(t) can be, at least in principle, 
experimentally detected are in general very small. A comparison with the exponential 
behaviour in (22), which gives 



V(t)~l-t/T = l-tTr[r Pw (0)] 



(38) 



indicates that this could happen only for times smaller than t m ~Tq/t = 1/(AE) 2 t. How- 
ever, by taking into account the Heisenberg time-energy uncertainty principle, AE At > 1, 
one finds t m to be actually smaller than the time interval At necessary to complete any 
measurement. [1, 23] Therefore, it is practically impossible to detect deviations from the 
standard exponential decay law at small times in any actual experiment involving elemen- 
tary particles (nevertheless, in suitable atomic systems the situation might be different, 
see [24, 25]). 

In studying the large time behaviour of U(t), it is convenient to rewrite (35) as: 



U(t) = [1 - R(0)] e~ iHwt + R(t) 



(39) 
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where 

1 _,YW(°) 



m = *<i> - ' < . ,40) 

For large times, the prefactor to the Wiesskopf-Wigner term, which is actually equal to 
dH\y(E )/dE , can be reabsorbed in a normalization of the states. By putting the system 
in a box and using again the definitions (25) and (26), in the infinite- volume limit the 
correction term R(t) takes the form 

poo -i(E-i\)t 

R (t) = J E JE. (E) (Eo _ E + iX)2 - (41) 

Due to phase space limitations, the behaviour on threshold of ui{E) can be usually approx- 
imated by: 

E ^ 



u(E)~u(E )[ — -l) , 5>0. (42) 



Inserting this expression in (41), together with a convergent factor e~ sE ', e <C t, the integral 
in R(t) can be evaluated exactly in terms of Whittaker functions. [26] Taking the large t 
limit, one finally obtains 



E m Lv(Eq) e -i(E m -i\)t 

(E -E m + t\) 2 (i t E m y 



ns + 1) -—stt (43) 



Therefore, for large times the effective evolution operator U(t) exhibits a power law be- 
haviour, modulated by an exponential. [22] Notice that the probability V(t) in (4) has a 
more complicated behaviour due to the interference effects between the two terms in (39) . 
In fact, besides the standard Wiesskopf-Wigner term proportional to e~*/ T and the power- 
like term proportional to t~ 2( - s+1 \ V(t) contains also an oscillating term that, taking for 
simplicity A = 0, is proportional to e~ t ' 2r He[e~ %tEm / {it) 5+1 \. In another context, this 
intermediate behaviour has also been noticed in [24]. 

In order to estimate the region in which the long-time correction (43) supersedes the 
standard exponential decay term, one has to consider the magnitude of the correction, 
given by the modulus \R(t)\ of the operator in (43). Recalling (296), one has 

m) ^r-2XAr\ ^ _ | :>2 - y ^ , ,44) 

where Ar = 1 — /j'(E ). This result should be compared with the standard exponential 
term e~ rt l 2 . In order to do this, let us choose a basis in which F (or better the spectral 
operator uj) is diagonal, and label with the index a the corresponding eigenvalues. In the 
case of the neutral kaon system, a takes the two values S and L, which refer to the K~s and 
Kl states. Assuming the difference E — E m of the same order of E m , one finds that the 
power law behaviour dominates for times larger than T a jr a , where r a is implicitly given 
by the following equation 

1 + 2 A) {S + 1) lnr a + 2 (l + 2 A) {S + 2) In f^j + 2A^ . (45) 
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Since E m j T a is usually very large (~ 10 15 ) in elementary particle decays, it turns out that 
deviations from the exponential law can be seen only after hundreds of life-times (for 6 = 1 
and A = 0), a region which is clearly unavailable to the experiment. Furthermore, notice 
that the role of the environment tends to worsen the situation, pushing the limit towards 
even longer times. 

We can safely conclude that the role played by the environment in the decay process 
of an unstable system is marginal and can not be detected in actual experiments. This 
conclusion results from a careful study of the starting evolution equation (7). Although 
phenomenological in nature, this equation encodes in a rather general and universal way 
possible dissipative effects due to a weak interaction with an environment. From this 
point of view, the description of the neutral kaon system in terms of completely positive 
dynamical semigroups given in [13-16] is appropriate. In particular, the deviations from 
the predictions of the standard Weisskopf-Wigner theory discussed there, if experimentally 
detected, could be really the sign of new physics. 
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